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Abstract 

The dependence of the multiplicity on the number of collisions and the trans- 
verse momentum distribution for central and peripheral Au-Au collisions are 
studied in the model of percolation of strings relative to the experimental 
conditions at RHIC. The comparison with the first RHIC data shows a good 
agreement. 
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The first experimental data from the relativistic heavy ion collider show several interest- 
ing features • The measured multiplicity lies in the lower range of the values predicted 
by various models and shows a smooth dependence on the number of participants [fj]||. The 
experimental pt distributions in Au-Au central collisions reveal a large departure from the 
naive perturbative quantum cromodynamics (pQCD) predictions Also, some differ- 

ences are seen for peripheral collisions. In this paper, we compare these experimental results 
with the predictions deduced from the percolation of strings. 

Multiparticle production is currently described in terms of color strings stretched be- 
tween the projectile and target, which decay into new strings and subsequently hadronize 
to produce observed hadrons. Color strings may be viewed as small areas in the transverse 
space, TirQ, r =0.2 fm, filled with color field created by the colliding partons PJT5I- Particles 
are produced via emision of qq pairs in this field. With growing energy and/or atomic num- 
ber of colliding particles, the number of strings grows, and they start to overlap, forming 
clusters, very much like disks in the two-dimensional percolation theory |9|{ll|. At a certain 
critical density a macroscopic cluster appears which marks the percolation phase transition. 
The influence of the clustering of strings and their percolation on the transverse momentum 
spectra was studied in in the thermodynamical limit of a very large interaction area S 
and very large number N of formed strings, corresponding to very high energies and atomic 
numbers of participants. In this limit the study allows for an analytical treatment. Here 
we apply this approach to realistic energies and colliding nuclei at RHIC using Monte- Carlo 
simulations. 

The percolation theory governs the geometrical pattern of the string clustering. Its 
observable implications, however, require introduction of some dynamics to describe string 
interaction, i.e., the behaviour of a cluster formed by several overlapping strings. 

There are several possibilities [|IT|-|I"3|j . Here we assume that a cluster behaves as a single 
string with a higher color field Q n corresponding to the vectorial sum of the color charge 
of each individual Qi string. The resulting color field covers the area S n of the cluster. 
As Q n = J2iQi, an d the individual string colours may be oriented in an arbitrary manner 
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respective to one another, the average QuQij is zero, and Q\ = nQ\. 

Knowing this charge color Q n , one can compute the particle spectra produced by a single 
color string of area S n using the Schwinger formula |L4||. For the multiplicity fi n and the 



average of particles, < p\ > n , produced by a cluster of n strings one finds 



/ n Sn o / 2 / 1 \ 

tin= \~ST fJ ' 1 ' < Pr>n= W— <P T >i I 1 ) 

where \i\ and < p\ >\ are the mean multiplicity and p\ of particles produced by a simple 
string with a transverse area S\ = tttq. For strings just touching each other S n = nS\, 
and hence \i n = n/xi, < p\ > n —< p\ >i, as expected. In the opposite case of maximum 
overlapping S n = Si and therefore \i n = s/nfii, < p\ > n = y/n < p\ >i> so that the 
multiplicity results maximally supressed and the mean p\ maximally enhanced. Notice that 
a certain concervation rule holds for the total < p\ >: \i n < p\ > n = < p\ >\. 

Eq. (HD is the main tool of our calculation. In order to compute the multiplicities we 
generate strings according to the quark-gluon string model and using the Monte-Carlo code 



described in ]TS]. Each string is produced at an identified impact parameter. From this, 
knowing the transverse area of each string, we identify all the clusters formed in each collision 
and subsequently compute for each of them its multiplicity in units The value of \i\ was 



fixed by normalizing our results to the SPS WA98 results for central Pb-Pb collisions |L6 
Note that in colour string models both \i\ and < p\ > are assumed to rise with energy due 
to the increase of the rapidity interval available and hard scattering contribution. However 
this rise is very weak in the the energy range \fs = 17 -7- 200 GeV. Therefore we can neglect 
the energy dependence of both fix and < p\ >% in the first approximation. 

The comparison of our results for the dependence of the multiplicity on the number of 
participants with the SPS WA98 data, with the RHIC PHENIX § and PHOBOS data 
at v/i = 130 GeV and with PHOBOS data at y/s = 200 GeV is presented in Fig. |. The 
agreement with the data at yfs = 130 GeV is quite good. Our predictions for yfs = 200 
GeV lie a little below the data. This could be a consequence of the mentioned weak energy 
dependence of p,\. On the whole the agreement with the data indicates that the percolation 



of strings mechanism correctly describes the behaviour of multiplicities in nucleus-nucleus 



collisions. This conclusion has also been reached in p0|. 

In order to compute the transverse momentum distributions, we make use of the 
parametrization of the pp UAl data at 130 GeV by A. Drees || 

— = — (2) 

dvl (po+p T ) a 

where a,po and a are parameters fitted to data [[J. The standardly expected distribution 
for central Au-Au collision is 



dN 
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However, if we take into account percolation of strings then from Eq. p] we conclude 



where, by definition, 
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Here the sums go over all clusters formed in Au-Au collisions and fi n are given by Eq. (1) 
To exclude the unknown < p\ > x we compare Eq. (4) with the one for pp collisions 



< Pt > pp =< ^ >l' p 2 < P 2 T >i (6) 

to finally get 
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This equation implies that the same parametrization @ for the transverse momentum dis- 
tribution can be used for nucleus-nucleus collisions, with the only change 

/ < nS\/S n >Au-Au\ ^ / c x 

Po^Pol - g /Q ^ (8) 

V < nSx/S n > pp J 

Note that at the energies considered only two strings are exchanged in pp collisions on the 
average, with very small fusion probability, so that < nS\/S n > pp — 1 



In Fig. we show the distribution (|2|) with the change (Q) for central (5%) Au-Au 
collisions (solid line) compared to the PHENIX experimental data P,|IT] (black squares). 
Also the distribution expected from the independent string picture (Eq. (3)) is shown 
(dash-dotted line). For peripheral (80-92%) Au-Au collisions the same results are presented 
in Fig. |3| (with the same notations). A very good agreement is observed in both cases. We 
stress that this result is obtained using a very simple formula deduced in a simple way from 
the dynamics, with practically a single parameter, the string transverse radius. 

Our predictions for central (5%) Ag-Ag and central Au-Au collisions at 200 GeV/n are 
presented in Fig. 0. 



In the thermodynamic limit one approximately obtains [12 



< ^ >= ; — = v*^) (9) 

S n 1 - exp(-?7) 

where rj = Nnr^/S is the percolation parameter assumed to be finite with both iV and S 
large. Then one obtains an analytic expression for the transverse momentum distribution 
for all colliding particles in terms of rj 

dN a 

dp 2 t ~ (p ^/F( Vpp )/F(ri)+ PT ) a 

In (10) a and po are parameters of the experimentally observed spetrum in pp collisions 
at a given energy and rj and r\ pp are percolation parameters for nucleus-nucleus and pp 
collisions respectively, to be obtained from the geometrical percolation picture. Naturally 
all parameters depend on energy, although, as mentioned, at energies below 200 GeV rj pp is 
quite small so that F(rj pp ) ~ 1. Eq. (10) implies a certain universality of the distributions 
at fixed energy, for all participant nuclei and different centralities: the distribution results a 
universal function of r). 

At very high energies both rj and r\ pp become large, so that F{r\) ~ 1/ \ff\. Then we find 

F( V ) ' ' 

For minimum bias collision of identical nuclei 




, = 16^ (12) 

a AA 

where v the average number of collisions and we used the fact that the average overlap 
area is 1/16 of the total nucleus- nucleus cross-section. From this we find an asymptotic 
expression for the distribution valid at very high energies (large 77's): 

dN a (13) 



dPt (2p ^Aa pp /af A +p T ) a 
Accordingly we find for these energies 



< Pt >aa= 4:<p 2 T > pp (14) 
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<A 

from which it evidently follows that at very high energies the ratio of the averages p\ in 
nucleus-nucleus and pp collisions should grow as the total pp cross-section. 
Notice that from (1) we have for a cluster 

< P T >n= 7; (15) 

b n Hi 

In the limit of very high density of strings, as observed from the Monte-Carlo calculations, 
practically all strings overlap into a single cluster, which occupies the whole interaction 
area. This is obviously an alternative description of parton saturation. In this limit from 
(15) we find a universal relation between the observed < p T >ab and the multiplicity per 
unit overlap area in the (nucleus A)- (nucleus B) collisions, valid for any participant nuclei 
and any centrality: 

/ 2 ^ Si <Pr >i 1 / 1R N 

< P T >ab= 7 — fi>AB (16) 

A*l ^AB 

Here Sab is the overlap area depending on the impact parameter. It is remarkable that a 
similar relation has been found using the idea of the saturation of gluons in [25]. In their 
case the proportionality coefficient is given by a s P 2 (m h /p s ), where m h is the mass of the 
produced hadron p s is essentially the saturation momentum and P is some function which 
tends to a constant as its argument vanishes. At high energies and therefore large p s this 
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coefficient tends to a constant modulo logarithmic terms coming from the running of the 
coupling constant. In our case the coefficient also tends essentially to a constant at high 
energies, since both Si < p\ > and fi\ (per unit rapidity) tend to a constant. 

To conclude, some general remarks. The colour string model, in principle, was introduced 
to describe the soft part of the spectra. However our results remain in good agreement with 
the experimental data up to pt of the order 5 GeV/c, where hard scattering is supposed to 
play the dominant role, standardly desribed by the PQCD approach. This implies that with 
the parametrization (||D the hard part of the spectrum is also included into our string picture, 
corresponding to the high-momentum tail of the string spectra. Clustering and percolation 
of strings can then be looked upon as an alternative description of non-linear effects, which 
in the PQCD language correspond to high partonic densities, such as saturation of partons. 

There are several models which explain qualitatively the dependence of the multiplicity 
on the number of participants [fL8[ - |23[| , although not all of them are able to reproduce quan- 
titatively the experimental data. The situation is worse concerning the pt distributions. 
With some fitted parameters, the jet quenching effect [§,[23J] could explain the data, but in 
this case one expects an enhancement of the multiplicity in the central rapidity region, which 
is not seen by experimental data. The percolation of strings is a natural and simple way of 
explaining both the multiplicity and transverse momentum distribution. To further check 
our approach, experimental information on forward-backward correlations of multiplicity 



and transverse momentum distributions would be welcome (see 11,121). 
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FIGURES 
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FIG. 1. Comparison of the dependence of the multiplicity on the number of participants with 
the SPS WA98 [16] data (filled triangles), with the RHIC PHENIX [5] (filled boxes) and PHOBOS 
[6] (non-filled boxes) data at yfe = 130 GeV and with PHIC PHENIX data at ^200 GeV (filled 
stars) . The dashed, solid and dotted lines are our predictions for the relevant energies. 
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FIG. 2. Expected px distribution using the percolation of string model (solid line) for central 
(5%) Au-Au collisions compared with PHENIX experimental data [2,17] (filled boxes). Also it is 
shown the expected distribution [8] given by formula (||) and (||) (dotted-dashed line). 
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FIG. 3. Expected px distribution using the percolation of string model (solid line) for peripheral 
(80-90%) Au-Au collisions compared with PHENIX experimental data [2,17] (filled boxes). Also 
it is shown the expected distribution [8] given by formula (||) and (0) (dotted-dashed line). 



dN/ (7TNdp, 



10" 
10 : 



[GeV/c] 



10' 
10' 
10" 
10" 



3 



4 5 

p t (GeV/c) 

FIG. 4. Predictions for central (5%) Au-Au (solid line) and central (5%) Ag-Ag (dashed line) 
collisions at 200 GeV/n. 
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